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1 . Introduction. 



Motivated by applications to reliability theory, the model we shall 
consider is a mixture of the form 



F T (t) 



e At dG(A) , 

* 

0 



t ^ 0 



( 1 . 1 ) 



where F (t) is the survival function of the mixture, F(t;A) = e is the 

survival function of the component distribution given that A = A for A a 
nonnegative random variable having cdf (cumulative distribution function) G. 

We can refer to such a mixture as a G-mixture of exponentials. In the relia- 
bility context, T represents the time to failure of an item drawn at random 
from a population of items with different exponential life distributions, mixed 
according to the distribution of A. 

We take the usual definition for a failure rate function, i.e. for the 
random variable T having density f^,(t) and survival function F^,(t) , the 
failure rate function is defined by 

f T (t) 

h (t) = 3 , for all t such that F (t) > 0. 

F T (t) 

The initial mixing distribution G may be modified to yield a related 
distribution at any time t, i.e. let G^_(A) = P(AsA | T>t) and define A(t) 
to be a random variable having cdf G^ . G^ is the residual mixing distribu- 
tion at time t. It represents the revised distribution for the failure rate 
of a randomly drawn item from the mixture, given that the item is still alive 
at time t. The relationship between G and G may be exhibited through a 
straightforward application of Bayes' Theorem which yields 



dG(A), 



( 1 . 2 ) 



-At 



dG (A) = ^ r 
E(e At ) 



where the expectation is with respect to G, the distribution of A. 

For such a mixture of exponentials as (1.1), it is well known (Barlow, 
Marshall, and Proschan [2]) that the mixture failure rate h^(t) is decreasing in 
(we shall use ’decreasing 1 for non-increasing and ’increasing’ for non-decreas- 
ing), and that (Gnedenko, Belyayev, and Solovyev [3]) the mixture failure rate 
function approaches as a limit with increasing time t the least of all 
parameter values positively present in the mixture, i.e. A^. 

Aldrich and Morton [1] and O’Bar [4] showed that at time t the failure 
rate function of a mixture of exponentials is the expected value of the resid- 
ual mixing distribution corresponding to that time, i.e. 



h T (t) = E [A(t) ] (1.3) 

Thus, the sequence of residual mixing distributions is decreasing in expectation. 
This fact supports the intuitive notion that with increasing time the shorter- 
lived members of the mixture are being "weeded out." 

The conjecture of a compatible stochastic ordering of the residual 
mixing distributions arises as a logical extension of the foregoing results. 

In this paper the existence of such a stochastic ordering is demonstrated. It 
is also shown that the residual mixing distributions converge in distribution 
to the distribution degenerate at the least parameter value positively present 
in the mixture, with an attendant monotone convergence of raw moments. The 
latter result is one asserted by Aldrich and Morton. 



2 



The results are presented as applications of Propositions 1 and 2, to 
be developed in the next section. These propositions state properties of 
E{u[A(t)]} for fairly general functions u. To indicate the relevance of 
these propositions, we first note that an expression of this form is involved 
in (1.3), taking u[A(t)] = A(t). We next define an indicator function 



! 0, ASX 

1, A > X 

Using this indicator' function, may be expressed as 



G t (X) 



dG fc (x) = EU x [A(t)]}, 
X 



another expression of the same form, with u[A(t)] = I^[A(t)]. 



(1.4) 



(1.5) 



2. Basic Results . 

We shall develop in Propositions 1 and 2 of this section basic tools 
which should prove to be useful in a wide variety of applications. Several 
of these applications are given in section 3. 

Using (1.2) we find that 



E{u[A(t) ] } 



u(X)dG fc (X) 



u(^)e 



-Xt 



E(e At ) 



dG(X) 



ETu(A)e At 1 
E(e“ At ) 



( 2 . 1 ) 
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Thus, expectations of functions of A(t) may be evaluated by simply taking 
expectations of functions of A. A straightforward application of the Monotone 
Convergence Theorem yields 



{E[u(A)e At ]} = E[u(A) ^ e' At ] 
= -E[Au(A)e" At ], 



( 2 . 2 ) 



for A £ 0. 

We next require the function u to be nonnegative and monotone to get 



PROPOSITION 1 . If A is a nonnegative random variable and u a nonnegative 
monotone increasing (decreasing) function, then 



-At-, 



<P„(t) = E[ - U(A ^r. ] = E{u[A(t) ] } 

E(e‘ At ) 



u 



is monotone decreasing (increasing) in t. 



Proof : Using (2.2), 



d . E(e~ At )E[-Au(A)e~ At 1 - Efu(A)e~ At lEr-Ae~ At 1 

dt " u(t) ' E ( e - At ) 2 



E [Au(A)e At 1 + E [u(A)e AT 1 E(Ae~ At ) 

-At. -At. -At. 

E(e ) E(e ) E(e ) 



= E{u[A(t)]}E[A(t)] - E{A(t)u[A(t)]} 



= - C0V{u[A(t) ] , A(t) } . 



Then by the proporty of the covariance of similarly/dissimilarly ordered 
functions of the same random variable. 
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(i) u increasing in A(t) =* — cp^Ct) ^ 0, and 



(ii) u decreasing in A(t) =» — cp^Ct) ^ 0. 



□ 



We now need two lemmas regarding limits of expectations involving 
time t as a parameter. 

LEMMA 1 . Given u a nonnegative function. 



lira E[I x (A)u(A)e (A X)t ] = 0, 



Proof : Note that I^CA) = 1 ** A - A > 0 « e ^ A X ^ C 4- in t. An application 



of the Monotone Convergence Theorem yields 



! 



lim E[I x (A)u(A)e (A X)t ] = E[I (A)u(A) lim e (A X)t ] =0. □ 

t "^ 00 



LEMMA 2. For A > X = inf{A | G(A)<1}, 

o 

lim E [ (1-1. (A))e” (A " X)t ] > 0. ' 

Proof: Note that (1-1 (A)) =1^A-A^0« e -( A_x )t ^ 1 , f or t ;> 0 . Then 

A 

E[(l-I. (A))e' (A_X)t ] JE[1-I, (A)] = G(A) , and lim E[(l-I. (A) )e _(A_X)t ] * G(A) > 0. 

□ 

If, in addition to the mono tonicity requirement of Proposition 1, we 
also require u to be right-continuous at A^, we may find the limiting value 
of (p (t) as t -> 00 . We state this result in 



PROPOSITION 2 . . Let A be a nonnegative random variable having cdf G, and 
u a nonnegative monotone increasing (decreasing) function. If u is right- - 
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continuous at X^ = inf{X | G(X)<1), then cp^(t) decreases (increases) 
monotonically to u(X^) as t -> °°. 



Proof : 

(i) If u(A) + in A then cp^(t) ^ u(X q ). Proposition 1, <p^(tH in t. 

Thus, lim cp (t) exists and lies in the interval [u(X ) ,°°) . Choose 
u o 

t-x» 

X > X . Then, 
o 

E[(l-I x (A))u(A)e _At ] + E[I x (A)u(A)e" At ] 

cp (t) = — — 

U E[(l-I x (A))e At ] + E[I x (A)e" At ] 

E[(l-I x (A))u(A)e' (A_X)t ] + E[I x (A)u(A)e _(A_X)t ] 
E[(l-I x (A))e' (A_X)t ] + E[I x (A)e' (A_X)t ] 

Since u(A) + in A, 



<p u (t> ^ 



u(A)E[(l-I x (A))e (A X)t ] + E[I x (A)u(A)e (A X)t ] 
Et(l-I x (A))e _(A_X)t ] + E[I x (A)e _(A_X)t ] 



^ u(X) + 



E[I x (A)u(A)e (A X)t ] 



E[(l-I x (A))e (A X)t ] + E[I x (A)e (A X)t ] 



( 1 ) 



Lemmas 1 and 2 may be applied to show that the last term in expression 
(1) goes to zero as t Therefore, lim cp (t) ^ u(X) . Let X I X 

t“Xx> u 

to get lim cp (t) £ u(X ). Then, cp (t) + u(X ) as t + 00 . 

t“H» u u 

(ii) The proof for u monotone decreasing is similar to that above. □ 



3 . Applications . 

The first application of the tools developed in the preceding section 
concerns the stochastic ordering of the sequence of residual mixing distribu- 
tions and its limiting distribution, stated as 
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Theorem 1 . If A is a nonnegative random variable with cdf G and 



X = inf{X I G(A)<1}, then for a G-mixture of exponentials 
o ’ 

a) G^ is a stochastically decreasing sequence in t, i.e. 

t- ^ t~ ^ G (X) ^ G (X) , for all X • 

1 1 C i t 2 

b) G fc converges in distribution to the distribution degenerate at 

X , as t -> °°. 
o 

Proof : 

a) By (1.5), 

G t (X) = E{I x [A(t)]}. 

Note that 1^ is monotone increasing in its argument. Application 

of Proposition 1 with u[A(t)] = I, [A(t)] yields G (X) + in t. 

A t 

b) Application of Proposition 2 with u[A(t)] = I^[A(t)] yields as 
t -*■ 00 , 

S O, X ^ X 

o 

1, X > X . □ 

o 

An immediate consequence of Proposition 2 is the well-known result that 
a mixture of exponentials is DFR, with its failure rate approaching the least 
parameter value positively present in the mixture. We state this result as 

J 

Theorem 2 . If A is a nonnegative random variable with cdf G and 

X q = inf{X | G(X)<1}, then for a G-mixture of exponentials the mixture failure 

rate function, h^,(t), decreases monotonically to X q as t 00 . 
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Proof : As noted in (1.3), 



h T (t) = E [ A ( t ) ] . 

Application of Proposition 2 with u[A(t)] = A(t) yields the desired result. 

□ 

We may also see the monotone convergence for the raw moments of the 
residual mixing distribution in 

Theorem 3 . For p ^ 0, E[A(t)^] decreases monotonically to X^ . 

Proof : Application of Proposition 2 with u[A(t)] = A(t)^ yields the desired 

result. □ 

The variance of A(t) is of interest, since it has an interpretation 
as the negative of the slope of the mixture failure rate function (O'Bar [4]). 
We state the following as 

Corollary 1 . The limit as t 00 of VAR[A(t)] is zero. 

Proof: We express the variance of A(t) as 

VAR[A (t) ] = E[A(t) 2 ] - E 2 [A (t) ] (3.1) 

Then, by Theorem 3, 

lira VAR[A(t) ] = X 2 - X 2 = 0. □ 

o o 

t-X» 

Expression (3.1) leads us to 

Corollary 2 . On the interval (0,«>), VAR[A(t)] is of bounded variation. 

Proof ; From the power series expansion for e^ , we have that 

,n n! At 
A < — e , 
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for n a nonnegative integer. Thus, 



,.n -At. , At -At. 
E [A(t) n ] - E i A « 1 < E < e « 1 

E[e C ) t E(e C ) 



< 



n! 



t n E(e“ At ) ' 



for t > 0. 



Then, from expression (3.1) and Theorem 3, since VAR[A(t)] is expressed as 
the difference of two real-valued nonnegative monotone functions for t > 0, 
it is of bounded variation on (0,°°). □ 
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